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Abstract 

A  clean  hydrogen  (99.9999%)  and  dry  ice  production  process  is  proposed,  which  is  based  on  phase  equilibrium  (PE)  separation  and 
methane  reforming.  Heat  and  power  integration  studies  are  carried  out  for  the  proposed  process,  by  formulating  and  solving  the  minimum 
hot/cold/electric  utility  cost  problem  for  the  associated  heat  exchange  network.  The  optimum  operating  cost  of  the  proposed  process  is  shown 
to  be  lower  than  the  corresponding  cost  of  the  conventional  PSA  (pressure  swing  adsorption)  based  process,  if  the  produced  dry  ice  is  sold 
for  as  low  as  2  cents  kg-dry-ice^1  or  if  an  equivalent  CO2  sequestration  credit  is  conceded. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  highly  efficient  oxidation  of  hydrogen  to  water  in  fuel 
cells  suggests  it  as  an  environmentally  attractive  transporta¬ 
tion  fuel,  whose  use  could  result  in  a  positive  health  impact 
on  city  populations  [1  ].  The  most  common  industrial  process 
for  production  of  hydrogen  from  natural  gas  is  steam  reform¬ 
ing  [2-4],  which  involves  the  endothermic  transformation  of 
methane  and  water  to  hydrogen,  carbon  dioxide  and  carbon 
monoxide.  Following  its  formation,  hydrogen  must  be  sep¬ 
arated  from  the  other  gases,  to  attain  purity  levels  required 
for  best  performance  and  long  operation  of  fuel  cells  [5]. 
This  is  especially  the  case  for  proton  exchange  membrane 
fuel  cells  (PEMFC),  whose  anode  platinum  catalyst  has  an 
extremely  low  carbon  monoxide  (CO)  tolerance  (only  few 
ppm)  [6].  Hydrogen  concentration  is  increased  prior  to  sepa¬ 
ration,  through  the  use  of  water  gas  shift  reactors  where  CO 
is  partially  consumed  down  to  a  low  percentage  as  it  reacts 
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with  water  to  produce  more  Ht  and  CCH.  Most  of  the  remain¬ 
ing  water  is  later  separated  by  condensation.  CO  content  can 
sometimes  be  reduced  to  few  ppm  either  by  using  methana- 
tion  [7]  reactors,  i.e.  catalyzing  the  reaction  of  CO  with  Ht  to 
produce  CH4,  or  through  CO  preferential  oxidation  (PROX) 
[8]  reactors,  which  require  the  addition  of  precisely  measured 
amounts  of  air. 

Conventional  hydrogen  production  processes,  however, 
[4,9-12]  use  pressure  swing  adsorption  (PSA)  technology 
for  final  hydrogen  purification  [13,14].  PSA  achieves  sepa¬ 
ration  of  CO,  CO2,  CH4  and  H2O  from  H2,  by  adsorption 
of  these  components  on  a  solid  adsorbent  at  a  relatively  high 
pressure.  The  adsorbed  species  are  then  desorbed  from  the 
solid,  by  lowering  the  pressure  and  purging  with  high  purity 
product  hydrogen.  The  resulting  PSA  waste  gas  contains  sig¬ 
nificant  amounts  of  hydrogen  and  methane  and  is  thus  burned 
as  a  source  of  heat  for  the  reformer.  Continuous  flow  of 
hydrogen  product  is  maintained  by  using  multiple  adsorp¬ 
tion  beds,  whose  adsorption/desorption  cycles  are  properly 
synchronized. 

In  this  work,  a  clean  hydrogen  (99.9999%)  and  dry  ice 
production  process  is  proposed,  which  is  based  on  phase 
equilibrium  (PE)  separation  and  methane  reforming.  The  pro- 
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Nomenclature 

Ccvcj 

cost  coefficient  of  cold  utility  with  constant 
temperature  number;  =  1,  2,  3  ($kJ-1) 

CCUVj 

cost  coefficient  of  cold  utility  with  varying 
temperature  number;  =  1,  2  ($kg-1) 

CHU 

cost  coefficient  of  hot  utility  ($  kg-1) 

Cp 

mass  heat  capacity  (kJkg-1  K-1) 

CW 

cost  coefficient  of  electric  utility  ($  kJ-  1 ) 

c 

dry  ice  selling  price  or  CO2  sequestration 
credit  ($  kg -dry-ice'  1 ) 

cue j 

cold  utility  with  constant  temperature  number 
7=  1,2,  3 

CUV; 

cold  utility  with  varying  temperature  number 
7=1,2 

Den 

denominator  in  reaction  rate  expressions 

F 

mass  flow  (kg  s-1) 

HE 

heat  exchanger  subnetwork 

HEP 

heat  engine  and  pump  subnetwork 

HU 

hot  utility 

AH 

enthalpy  change  (kJ  s-1) 

A 

standard  (25  °C,  1  atm)  heat  of  reaction  rm, 
m  =  1,  2,  3  (kJ  mol-1) 

kk 

adsorption  constant  of  species  k  =  CH4,  H2O, 
H2,  CO.  Units  are  specified  in  Table  A.l 

km 

rate  coefficient  of  reaction  rm,  m  =  1 , 2, 3 .  Units 
are  specified  in  Table  A.  1 

Km 

equilibrium  constant  of  reaction  rm,  m  =1,2, 

3.  Units  are  specified  in  Table  A. 2 

MUC 

minimum  utility  cost 

Pk 

partial  pressure  of  species  k  =  CH4,  H2O,  H2, 
C02,  CO  (bar) 

PE 

phase  equilibrium,  PE  based  process 

PSA 

pressure  swing  adsorption,  PSA  based  process 

Q 

heat  flow  (negative  for  cooling)  (kJ  s-1) 

AS 

entropy  change  (kJ  K- 1  s- 1 ) 

I'm 

reaction  m  =  1 ,  2,  3 

>'rm 

rate  of  reaction  r,„,  m=  1,  2,  3 

(kmolkgcat-1  h-1) 

SMR 

steam  methane  reformer 

T 

temperature  (K) 

rr  H 

1  i 

high  temperature  of  interval  i  in  the  hot  tem¬ 
perature  scale  (K) 

2i+ 1 

low  temperature  of  interval  i  in  the  hot  temper¬ 
ature  scale  (K) 

A  Train 

minimum  approach  temperature  (K) 

w 

electricity  (negative  if  produced)  or  work  (neg¬ 
ative  if  work  is  done  by  the  fluid)  (kJ  s-1) 

ws 

work  provided  to  the  HEP  subnetwork  (nega¬ 
tive  if  work  is  produced)  (kJ  s- 1 ) 

Greek  letters 

cti 

variable  indicative  of  presence  of  hot  utility  in 
interval  i 

Yji 

variable  indicative  of  presence  of  cold  utility 
with  varying  temperature  number  j  in  interval 

ij=  1,2 

St 

available  heat  at  interval  i  (kJ  s- 1 ) 

SaJCj 

heat  transferred  to  cold  utility  CUC;  (kJ  s-1) 

m 

fraction  of  cold  composite  stream  from  interval 
i  used  in  the  HE  subnetwork 

0t 

fraction  of  hot  composite  stream  from  interval 
i  used  in  the  HE  subnetwork 

variable  indicative  of  presence  of  cold  utility 
with  constant  temperature  number ;  in  interval 
«',7'=  1,2,  3 

Subscripts 

C 

cold  composite  stream 

cat 

catalyst 

CUC; 

cold  utility  with  constant  temperature  number 
7=  1,2,  3 

CUV; 

cold  utility  with  varying  temperature  number 

7=  1,2 

H 

hot  composite  stream 

HU 

hot  utility 

i 

interval  i  for  optimization  problem 

in 

inlet  temperature 

j 

utility  number  for  optimization  problem 

n 

number  of  intervals  for  optimization  problem 

out 

outlet  temperature 

W 

electric  utility 

Superscripts 

C 

cold  temperature  scale 

cue ] 

cold  utility  with  constant  temperature  number 

7=  1,2,  3 

CUV; 

cold  utility  with  varying  temperature  number 

7=1,2 

H 

hot  temperature  scale 

HU 

hot  utility 

posed  flow  diagram  avoids  the  following  limitations  of  the 
PSA  based  process:  (1)  close  to  12%  [12]  of  the  hydrogen 
produced  in  the  reformer  and  the  water  gas  shift  reactors  ends 
up  in  the  PSA  waste  gas.  (2)  A  considerable  percentage  of  the 
methane  fed  into  the  reformer  is  not  converted  to  hydrogen 
resulting  in  increased  flows  through  all  reaction  and  separa¬ 
tion  units;  such  percentage  depends  on  the  reformer  operating 
conditions  and  it  is  around  20%  [  1 2]  for  pressure  of  25 .7  atm, 
outlet  temperature  of  1 130  K  and  steam/CPL;  molar  ratio  of 
3.12  in  the  feed.  (3)  The  unsteady  operation  of  the  adsorp¬ 
tion  (PSA)  beds  necessitates  the  use  of  precise  and  complex 
control  systems  to  ensure  continuous  flow  of  clean  hydro¬ 
gen.  In  addition,  heat  and  power  integration  studies  of  the 
proposed  PE  based  process  are  carried  out,  with  the  purpose 
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of  finding  the  minimum  hot/cold/electric  utility  cost.  Use  of 
the  same  methodology  previously  applied  to  the  conventional 
PSA  based  process  [12],  allows  comparison  of  operating  cost 
and  CO2  emissions  for  both  processes  and  leads  to  our  final 
conclusions. 


2.  Process  description 

The  proposed  hydrogen  and  dry  ice  production  process, 
based  on  phase  equilibrium  (PE)  separation  and  methane 
reforming,  can  be  represented  with  the  block  diagram  in 
Fig.  1.  Water  and  methane  at  ambient  conditions  are  fed  into 
this  process  and  transformed  into  clean  hydrogen  (99.9999%) 
and  dry  ice  (98.5%).  Chemical  transformations  only  take 
place  in  the  steam  methane  reformer  (SMR)  and  are  rep¬ 
resented  by  the  reversible  reactions  below: 

CH4  +  H20  =  CO  +  3H2  A  Hi  :  206. 1  kJ  mol-1  (n) 

CO  +  H20  =  C02  +  H2  AH2°  :  -41.15  kJmoU1  (r2) 

CH4  +  2H20  =  C02  +  4H2  A Hj  :  164.9  kJ  mol-1  (r3) 

The  kinetics  of  these  reactions  on  a  Ni/MgAl204  catalyst 
have  been  studied  by  Xu  and  Froment  [15]  and  are  listed  in 
Appendix  A.  Addition  of  reactions  ri  and  r2  results  in  reaction 
r3,  which  is  the  overall  steam  methane  reforming  reaction, 
and  thus  the  overall  reactor  operation  requires  that  heat  be 
provided  to  the  reformer.  Hydrogen  is  obtained  together  with 
all  the  other  species  and  a  purification  process  is  therefore 
required.  A  series  of  phase  equilibrium  based  separation  units 
are  used  here  to  sequentially  separate  H20,  C02,  CH4  and 
CO,  and  to  obtain  a  clean  hydrogen  gas  stream.  Water  is 
first  separated  by  condensation  as  the  gas  stream  that  exits 
the  reformer  is  cooled  down  to  almost  ambient  temperatures. 
Further  cooling  to  163  K  allows  for  separation  of  C02  and 
generation  of  dry  ice.  Most  of  the  methane  is  removed  in  the 
CH4  removal  unit  A  at  77  K  and  recycled  into  the  reformer, 
after  appropriate  heat  exchange.  CO  is  absorbed  in  clean  liq¬ 
uid  methane  [16]  (process  raw  material)  that  is  fed  at  86  K 


into  an  absorption  tower  in  the  CO  removal  unit.  Coming  out 
of  this  unit,  in  gas  phase  at  170  K,  methane  and  the  recovered 
CO  are  sent  to  the  reformer.  CH4  and  CO  traces  remaining 
in  the  hydrogen-rich  stream  are  removed  by  condensation  at 
45  K  and  finally  clean  gas  hydrogen  with  99.9999%  molar 
concentration  is  obtained.  Heat  exchange  and  compression 
operations  are  used  to  deliver  hydrogen  at  ambient  temper¬ 
ature  and  300  atm,  which  is  the  typical  storage  pressure  for 
hydrogen  powered  vehicles  and  related  hydrogen-fueling  sta¬ 
tions. 

Different  than  the  conventional  methane  reforming  based 
hydrogen  production  process  [4,9-12],  our  proposed  flow 
diagram  does  not  use  a  PSA  system  for  hydrogen  purifica¬ 
tion  and  does  not  possess  any  of  the  conventional  process’ 
limitations  described  in  the  previous  section. 

Water  gas  shift  reactors,  typically  used  to  increase  hydro¬ 
gen  concentration  and  reduce  CO  content  of  the  gas  obtained 
from  the  reformer,  were  not  considered  here  since  our  only 
purpose  is  to  demonstrate  the  feasibility  of  the  proposed 
hydrogen  purification  process  for  a  flowsheet  that  only  con¬ 
tains  a  single  methane-reforming  unit.  Nevertheless,  similar 
gas  shift  reactor  containing  flowsheets  can  be  developed  and 
will  be  the  subject  of  future  research. 


3.  Process  simulation 

The  proposed  process  has  been  simulated  using  Aspen- 
Tech’s  process  engineering  software  HYSYS®  version  3.1 
[17].  Peng  Robinson’s  equation  of  state  was  used  as  the 
thermodynamic  fluid  package  for  this  simulation,  based 
on  AspenTech’s  recommendation  of  this  package  for  high 
hydrogen  content  systems  and  after  ensuring  that  the  process 
conditions  are  within  the  package’s  temperature  and  pressure 
applicability  ranges  [18]. 

The  HYSYS  process  flow  diagram  is  presented  in 
Figs.  1  and  2.  The  reformer  feed  is  set  to  25.7  atm  and 
81 1  K,  corresponding  to  values  within  typical  entrance  con¬ 
ditions  for  the  reformer  [9-11,19-21],  and  has  a  steam/CH4 


Fig.  1.  Process  flow  diagram  for  steam  methane  reforming  based  hydrogen  production  with  phase  equilibrium  (PE)  based  separation.  HYSYS  simulation. 
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Fig.  2.  Process  flow  diagram  for  mass  and  energy  flows  of  the  phase  equilibrium  (PE)  based  hydrogen  plant  simulated  in  HYSYS.  Heat  ( Q )  flows  are  shown 
as  positive  (+)  for  heating  and  (— )  for  cooling.  Work  ( W)  done  on  the  fluid  is  shown  as  (+).  Units  of  energy  are  k.l  s  1 . 


molar  ratio  of  3.12.  Excess  steam  is  used  to  reduce  by¬ 
product  carbon  formation  [4,11,19].  The  reformer  is  simu¬ 
lated  using  the  kinetic  models  proposed  by  Xu  and  Froment 
[15]  (see  Appendix  A).  Heat,  in  the  amount  of  288.89  kJ  s~', 
is  provided  to  the  reformer  to  maintain  a  high  reaction 
rate.  The  reformer  is  modeled  in  HYSYS®  through  a  series 
of  CSTRs  (continuous  stirred  tank  reactors):  twenty  2.36- 
kgcat  reactors  were  considered  in  this  study,  with  increas¬ 
ing  operation  temperatures  that  vary  from  825  to  1010  K. 
The  exit  mole  fraction  of  hydrogen  is  calculated  to  be 
0.34. 

Hydrogen  purification  starts  with  the  condensation  and 
flash  separation  of  liquid  water  from  the  reformer  exit,  by 
cooling  the  gases  to  almost  ambient  temperature  through  the 
removal  of  543.8  kJ  s_1 .  The  condensed  water  is  recycled  to 
the  reformer,  while  the  hydrogen  containing  gas  is  further 
cooled  to  163  K  (through  removal  of  51.62kJs_1)  allowing 
for  separation  of  CO2  and  generation  of  163.47  kg  h-1  of 
dry  ice  (98.5%  CO2).  The  resulting  hydrogen-rich  gaseous 
stream  is  then  further  cooled  to  77  K  (through  removal  of 
30.99  kJ  s-1 )  where  it  is  flash  separated  into  a  gaseous  stream 
with  98%  H2  and  a  liquid  stream  with  64%  methane.  The  last 
one  is  recycled  into  the  reformer  after  addition  of  72.2 1  kJ  s~ 1 
of  heat. 

Liquid  methane  in  the  amount  of  59.68  kg  h-1  is  gen¬ 
erated  in  a  liquefaction  unit  that  requires  70.50  kJs-1  of 
compression  work.  CO  is  absorbed  in  the  liquid  methane 
in  a  10-equilibrium  stage  absorption  tower  used  in  the  CO 
removal  unit,  reducing  CO  content  of  the  hydrogen-rich  gas 
stream  to  1  ppm.  Net  heat  of  11.12  kJ  s— 1  is  used  in  the  CO 
removal  unit  for  regeneration  of  the  liquid  methane  absorbent . 
Final  cooling  of  5.42  kJ  s_1  reduces  the  temperature  of  the 
gas  to  45  K,  low  enough  to  further  remove  CH4  and  CO 
by  condensation  and  obtain  29.67  kg  h-1  of  gas  hydrogen 
with  99.9999%  molar  concentration  (0.37  ppm  of  CH4  and 
0.12  ppm  of  CO).  The  clean  hydrogen  is  then  warmed  up  to 
ambient  temperature,  acting  as  a  cold  stream  for  heat  inte¬ 
gration,  and  pressurized  to  300  atm  with  the  consumption 


of  87.99  kJ  s  1  of  compression  work,  considering  85%  adia¬ 
batic  efficiency  compressors. 

4.  Heat  and  power  integration 

The  approach  proposed  by  Holiastos  and  Manousiou¬ 
thakis  [22]  for  the  calculation  of  the  minimum 
hot/cold/electric  utility  cost  for  heat  exchange  networks 
is  employed  for  the  heat/power  integration  of  the  above 
hydrogen/dry  ice  production  process.  The  problem  statement 
is  [22]:  given  a  set  of  process  streams  with  specified  flow 
rates,  inlet  temperatures,  and  fixed  outlet  target  temperatures; 
hot  and  cold  utility  streams  with  known  temperatures  and 
unit  costs;  and  electrical  (work)  utility  with  known  unit  cost; 
identify,  among  all  possible  heat  exchange/pump/engine 
networks,  the  minimum  total  (hot/cold/electric)  utility  cost 
necessary  to  accomplish  the  desired  thermal  tasks. 

For  the  application  at  hand,  the  set  of  process  streams 
is  defined  by  the  material  streams  in  the  process  flow  dia¬ 
gram  of  Fig.  2;  the  specifications  of  the  available  hot,  cold 
and  electric  utilities  are  defined  in  Table  1.  Methane  com¬ 
bustion  gas  is  considered  as  the  hot  utility  available,  whose 
inlet  temperature,  7®u  =  2168  K,  is  calculated  as  its  adia¬ 
batic  flame  temperature  for  110%  combustion  air  and  outlet 
temperature,  7)^  =  313  K,  is  estimated  for  emission.  Cool¬ 
ing  water  with  inlet  temperature  T®uvl  =  298  K  is  the  cold 
utility,  which  is  allowed  a  10°  increase  prior  to  emission,  fol¬ 
lowing  heat  exchange.  The  same  magnitude  of  temperature 
change  is  considered  for  cold  water  with  inlet  temperature 
yUUV2  _  278  K.  Cold  utilities  with  constant  temperature  are 
available  at  7™® 1  =  253,  7<:uc2  =  223  K,  7<:uc3  =  77  K  (liq¬ 
uid  nitrogen).  The  formulated  problem  therefore  includes 
both  utilities  with  constant  temperature,  as  Holiastos  and 
Manousiouthakis’  [22]  formulation  does,  and  utilities  with 
varying  temperature,  similar  to  Posada  and  Manousiouthakis’ 
[12]  formulation.  The  heat  and  power  integration  formulation 
for  the  case  where  hot  and  cold  utilities  are  not  allowed  to  be 
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Table  1 

Utilities  specification 


Utility2 

Tin  (K) 

Tout  (K) 

Cp.avgb  (kJkg-1  K-1) 

Cost  ($  kg  *) 

Cost  (SIlT1) 

HU:  CH4  combustion  gas 

7™  =  2168 

*5?  =  313 

cphu=  1-44 

chu  =  1.724  x  10~2c 

6.45  x  10”6 

CUV  1 :  cooling  water 

yCUVi _ 298 

Tg™  =  308 

Cpcuvi  =4.31 

ccuvi  =  8  x  10~5  [27] 

1.85  x  10"6 

CUV2:  cold  water 

rCUV2  _  278 

7-“V2  =  288 

cPcuv2  =  4.32 

ccuv2  =  8.6  x  10~4  [27] 

2  x  10~5 

CUC1 

jcuct  -253 

jcuct  =253 

- 

- 

ccuct=3.2x  10-5  [27] 

CUC2 

TcUC2  _  223 

TCUC2  _  223 

- 

- 

CCUC2  =6  X  Hr5  [27] 

CUC3:  liquid  nitrogen 

yCUC3  _  77 

yCUC3  _  77 

- 

- 

CCUC3  =  1-5X  10”4d 

W :  electricity 

- 

- 

- 

- 

cw=  1.25  x  10"5  [28] 

a  HU,  hot  utility;  CUV,  cold  utility  with  varying  temperature;  CUC,  cold  utility  with  constant  temperature;  W,  electric  utility  or  work. 
b  Average  mass  heat  capacity. 

c  Calculated  from  cost  of  0.3422  $  kg-1  for  CH4.  Cost  of  CH4  is  estimated  to  be  energetically  equivalent  to  cost  of  natural  gas  (0.3299  $  kg-1  [28]). 
d  Calculated  as  cost  of  electricity  required  to  liquefy  the  nitrogen  for  recycle  when  using  the  Stirling  gas  liquifiers  (SLG-1  or  SLG-4)  from  Stirling  Cryogenics 
&  Refrigeration  BY  [29]. 


used  for  work  generation  is  given  below: 


min  <;  chu^hu 

,0j,8i,  Fhu  ,  ffcuv  j ,  8cvc  j 


+  y>uv  jFcuvj+5~'ccuc  jScucj 

7=1  7=1 

+  CW  (  J2(Fcp)ci(Tf  -  7&l)(l  -  di) 

V  ;=  1 

-  -  ^iX1  -  w))  J  CD 


subject  to: 


Si  +  [(Fcp)ffir/i  +  /rHUtpHUa(](7(H  ~  Ti+0  -  si+l 


J^ScucjXj 

7=1 


(FcP)cA  +  J^FcuvjCPcuVjyji 

7=1 


X  (7f  -  j)  =  0  Vie[l,n] 


(2) 


Y.(Fcp)h/]  ~ 

=°  (3) 

Si  >0  V;  e  [2,  n]  (4) 

<5i  =  Sn+ 1  =  0  (5) 

0  <ru<\  Vie[l,n]  (6) 

0  <6i  <]  Vi  e  [1,  n]  (7) 


0  if  t?  < 

1  if  t"  <  r™ 


Vi  e  [1,  n] 


(8) 


J  0  if  T^+l  >  T, 

Yjl  ~  l  '  if  Tf+i  <  t: 

Vi  e  [1,  u],  V/  =  1,2 


CUV./ 

3Ut 

CUV  j 


or  7; 
and 


<  tcuvi 

—  ^  in 

yC  >  y  CUV  j 


(9) 


0  if  7^,  ^  7CUC( 
1  if  Tf+X  =  Tcvcj 


Vi  e  [1,  n],  Vy  =  1,2,3 

(10) 


Thu.  T’cuvi  ,  Tcuv2.  <5cuci ,  <5cuc2,  <5cuc3  >0  (11) 

Nomenclature  is  provided  in  a  separate  section  of  the 
paper.  The  objective  function  (1),  the  utility  cost,  is  the  sum  of 
the  costs  of  hot,  cold  and  electric  utilities.  Eq.  (2)  represents 
the  energy  balance  in  the  HE  (heat  exchanger)  subnetwork  for 
each  temperature  interval.  The  HE  subnetwork  is  built  with 
heat  exchangers;  fractions  rii  and  0j  of  the  enthalpies  from  the 
hot  and  cold  composite  curves  in  interval  i,  respectively,  are 
used  in  this  subnetwork.  The  respective  complementary  frac¬ 
tions  (1  —  r]i )  and  (1  —  0j)  are  used  in  the  HEP  (heat  engine 
and  pump)  subnetwork,  built  with  reversible  heat  engines  and 
heat  pumps.  Eq.  (3)  is  the  overall  entropy  balance  for  the  HEP 
subnetwork.  The  second  law  of  thermodynamics  is  expressed 
in  the  HE  subnetwork  by  the  fact  that  the  available  heat  at 
each  interval.  Si,  should  be  non-negative  (Eq.  (4)).  The  overall 
enthalpy  balance  in  the  HE  subnetwork  is  ensured  by  con¬ 
straint  5.  Solution  to  the  formulated  problem  provides  results 
for  heat  and  power  integration  of  the  process.  Heat  integra¬ 
tion  only  (no  power)  or  pinch  analysis  [23]  can  be  done  by 
considering  only  the  HE  subnetwork,  i.e.  by  setting  =  1  and 
&i  =  1  in  constraints  6  and  7,  respectively. 
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Table  2 


Heat  and  power  integration  results 


Utility 

Conven.  process 
PSAa  [12] 

Heat  integ.  PSA 
[12] 

Heat  and  power 

PSA  [12] 

Heat  integ.  PEb 

Heat  and 
power  PE 

HU(kJgH2-1) 

10.7 

0 

0 

57.8 

1.0 

CUV1  (kJgH2“‘) 

22.3 

24.2 

15.1 

42.4 

14.0 

CUV2(kJgH2-‘) 

0 

0 

0 

2.9 

0 

CUC1  (kJ  gH2_1) 

0 

0 

0 

0.4 

0 

CUC2  (kJgH2-‘) 

0 

0 

0 

0.8 

0 

CUC3  (kJgH2-‘) 

0 

0 

0 

2.2 

0 

tU-AkJgHj^1) 

6.2 

6.2 

-3.0 

19.3 

41.4 

MUC  ($kgH2 1)  min.  utility  cost 

0.19 

0.12 

-0.01 

1.14 

0.55 

a  PSA,  pressure  swing  adsorption  based  process. 
b  PE,  phase  equilibrium  based  process. 
c  W,  work  or  electricity  consumed  (negative  if  produced). 


dotted  line  represents  the  cold  composite  curve.  The  smallest 
temperature  difference  that  two  streams  leaving  or  entering 
a  heat  exchanger  can  have,  termed  the  minimum  approach 
temperature,  is  Armn  =  10  K.  It  is  attained  at  five  so-called 
pinch  temperatures:  453,  302,  288,  258  and  227  K.  The  use 
of  cooling  water  (CUV  1 )  in  the  amount  of  349.7 1  kJ  s~ 1  rep¬ 
resents  almost  the  whole  of  the  enthalpy  change  of  the  cold 
composite  curve  between  298  and  308  K  (it  looks  approx¬ 
imately  horizontal  in  Fig.  3  due  to  the  scale  of  the  plot). 
The  remaining  small  enthalpy  change  between  these  temper¬ 
atures  represents  water  and  recycle  stream  heating  prior  to 
their  entry  in  the  reformer.  The  other  cold  utilities  represent 
smaller  amounts  of  energy  whose  magnitudes  are  shown  in 
the  fourth  column  of  Table  2.  Hot  utility,  in  the  amount  of 
476.10kJs_1  or  57.8kJgH2_1,  is  required  here  to  provide 
heat  for  the  reformer,  different  than  the  PSA  based  process 
(see  Table  2)  where  sufficient  [12]  heat  is  obtained  after  com¬ 
bustion  of  the  PSA  waste  gas. 

The  enthalpy  diagram  of  the  optimal  HE  subnetwork  of  the 
PE  based  process  after  heat  and  power  integration  is  presented 
in  Fig.  4.  This  HE  subnetwork  does  not  use  cold  utilities  other 
than  cooling  water  at  298  K  (CUV1),  and  requires  only  a  very 
small  amount  of  hot  utility. 


2400 


2000 

1600 

T  (K) 

1200 

800 

400 


0 

0  20  40  60  80  100  120  140 

zlH,  Enthalpy  Change  (kj/s) 

Fig.  4.  Enthalpy  diagram  of  the  optimal  HE  subnetwork  of  the  methane 
reforming  based  hydrogen  production  process  with  phase  equilibrium  (PE) 
based  separation,  after  heat  and  power  integration.  AT^in  =  10  K. 


- Hot  Streams 

.  Cold  Streams 

AHHU  =  8.08  kj/s 


-dHCw,=  1 14.94  kJ/s 


5.  Results  and  discussion 

The  optimization  problem  is  solved  using  the  linear  pro¬ 
gramming  software  MINOS  5.5  [24].  The  calculated  min¬ 
imum  utility  cost  (MUC),  expressed  per  kg  of  hydrogen 
produced,  is  included  in  the  bottom  row  of  the  last  two 
columns  of  Table  2,  where  complete  results  are  summarized. 
The  first  three  columns  of  Table  2,  included  for  comparison, 
correspond  to  the  results  from  our  previous  study  [12]  on  the 
conventional  methane  reforming  based  hydrogen  production 
process  with  PSA  based  purification.  Specifically,  the  values 
in  the  first  column  of  Table  2  correspond  to  the  conventional 
PSA  based  process  without  any  integration  analysis  (no  opti¬ 
mization).  The  second  column  corresponds  to  the  results  after 
heat  (no  power)  integration  of  the  PSA  based  process  and  the 
third  column  corresponds  to  the  results  after  heat  and  power 
integration  of  the  same  process.  The  fourth  column  contains 
the  results  for  the  heat  (no  power)  integrated  PE  based  pro¬ 
cess  and  the  fifth  column  includes  the  results  for  the  heat  and 
power  integrated  PE  based  process. 

The  enthalpy  diagram  of  the  optimal  HE  network  of  the 
PE  based  process  after  heat  integration  is  presented  in  Fig.  3. 
The  solid  line  corresponds  to  the  hot  composite  curve  and  the 


AH,  Enthalpy  Change  (kj/s) 

Fig.  3.  Enthalpy  diagram  of  the  optimal  HE  network  of  the  methane  reform¬ 
ing  based  hydrogen  production  process  with  phase  equilibrium  (PE)  based 
separation,  after  heat  integration.  ATnim  =  10K. 
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Fig.  5.  Enthalpy  diagram  of  the  optimal  HEP  subnetwork  of  the  methane 
reforming  based  hydrogen  production  process  with  phase  equilibrium  (PE) 
based  separation,  after  heat  and  power  integration.  A7)njnc  =  A7irljnjf  =  5  K. 
Work  required  (Ws)  is  181.68  kJs-1. 

The  difference  in  the  enthalpy  change  of  the  cold  com¬ 
posite  curve  with  respect  to  that  of  the  hot  composite  curve 
in  Fig.  5  is  the  net  work  required  by  the  HEP  subnetwork 
(Ws  =  181.68  kJ  s-1).  Addition  of  this  amount  to  the  process 
work  requirements  (159.31  kJ  s-1  or  19.3  kJ  gH2-1)  gives  the 
total  work  needed  by  the  heat  and  power  integrated  PE  based 
process  (Table  2).  The  optimal  HEP  subnetwork  is  composed 
of  both  heat  engines  and  heat  pumps,  whose  correspond¬ 
ing  sections  can  be  visualized  in  the  entropy  diagram  of 
Fig.  6.  Those  sections  where  the  hot  curve  is  above  the  cold 
require  a  heat  engine,  while  those  where  the  cold  is  above  the 
hot  require  a  heat  pump.  It  is  easily  visualized  that  a  large 
need  for  heat  pumping  exists  on  the  optimal  HEP  subnet¬ 
work. 

After  comparison  of  the  results  shown  in  Table  2,  it  is  clear 
that  each  of  the  two  processes  (PSA  based  and  PE  based) 
find  their  lowest  utility  cost  after  heat  and  power  integra¬ 
tion  and  that  the  PE  based  process  has  a  much  higher  utility 
cost  due  to  its  operation  at  cryogenic  temperatures  to  pro¬ 
duce  clean  hydrogen  and  dry  ice.  The  economic  effects  of 
the  production  of  dry  ice  (for  sale  or  CO2  sequestration)  and 
the  better  use  of  raw  materials  attained  in  the  PE  based  pro¬ 
cess  is  taken  into  account,  in  addition  to  the  minimum  utility 


Entropy  Change  (kJ/K-s) 

Fig.  6.  Entropy  diagram  of  the  optimal  HEP  subnetwork  of  the  methane 
reforming  based  hydrogen  production  process  with  phase  equilibrium  (PE) 
based  separation,  after  heat  and  power  integration.  A7)llm(-  =  ATmm  n  =  5  K. 
Entropy  changes  are  balanced:  A  Ac  =  —  ASh- 

cost,  with  the  calculation  of  the  operating  cost,  as  shown 
in  Table  3.  The  PE  based  process  uses  raw  materials  water 
and  methane  in  approximately  the  stoichiometric  amounts 
required  by  the  overall  reaction  H3),  which  represents  reduc¬ 
tions  of  29%  in  process  methane  and  6%  in  process  water  with 
respect  to  the  PSA  based  process.  The  operating  cost  of  the 
PE  based  process  is  lowered  as  a  market  for  the  produced  dry 
ice  is  developed  or  a  CO2  sequestration  credit  is  conceded. 
Higher  dry  ice  price/credit  (C)  obviously  results  in  lower 
operating  cost  as  shown  in  Table  3  for  values  of  C=  16  and 
C  =  6  cents  kg-dry-ice- 1 .  It’s  thus  found  that  if  the  produced 
dry  ice  is  sold  for  as  low  as  6  cents  kg-dry-ice-1  or  equivalent 
CO2  sequestration  credit  is  conceded,  the  PE  based  process 
results  in  lower  operating  cost  than  the  heat  and  power  inte¬ 
grated  PSA  based  process.  Also  shown  in  Table  3,  a  value 
of  C  =  2  cents  kg-dry-ice- 1  results  in  an  operating  cost  lower 
than  the  conventional  (without  any  integration  analysis)  PSA 
based  process. 

Permanent  storage  of  carbon  dioxide  in  the  marine  envi¬ 
ronment  is  described  by  Murray  et  al.  [25],  who  propose 
a  technique  for  the  disposal  of  solid  carbon  dioxide  in  the 
form  of  dry  ice  blocks  shaped  as  torpedos  that  are  left  to 
fall  through  the  ocean  water  to  finally  penetrate  in  selected 


Table  3 

Operating  cost 


Resource 

Conven.  process  PSA 

Heat  integ.  PSA 

Heat  and  power  PSA 

Heat  integ.  PE 

Heat  and  power  PE 

Process  CH4  (kgCH4  kgH2-1) 

2.84 

2.84 

2.84 

2.01 

2.01 

Process  H2O  (kg^Okg^-1) 

4.83 

4.83 

4.83 

4.50 

4.50 

Dry  ice  product  (kgC02  kgH2-1) 

0 

0 

0 

5.51 

5.51 

Operating  cost2  ($kgH2-1) 

1.17 

1.10 

0.97 

1.83  —  Cx  5.51 

1.24- Cx  5.51 

If  C  =  0. 1 6b  ($  kg-dry-ice- 1 ) 

1.17 

1.10 

0.97 

0.95 

0.36 

If  C  =  0.06b  ($  kg-dry-ice-1) 

1.17 

1.10 

0.97 

1.50 

0.91 

If  C=0.02b  ($  kg-dry-ice-1) 

1.17 

1.10 

0.97 

1.72 

1.13 

a  Operating  cost  =  MUC  +  process  CH4  cost  +  process  H2O  cost  — dry  ice  credit.  Process  water  price  of  5.3  x  10-4$kg-1  [27]  and  methane  price  of 
0.3422  $  kg-1  were  used. 

b  C:  dry  ice  selling  price  or  CO2  sequestration  credit  ($  kg-dry-ice-1). 
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Table  4 

CO2  emissions 


Source 

Conven.  process  PSA 

Heat  integ.  PSA 

Heat  and  power  PSA 

Heat  integ.  PE 

Heat  and  power  PE 

Plant  CHp  (kgCH4  kgH2~1) 

3.05 

2.84 

2.84 

3.10 

2.03 

Plant  C02b  (kgC02kgH2-‘) 

8.25 

7.71 

7.71 

2.99 

0.05 

Work  associated  C02c  (kgC02  kgHq1) 

0.74 

0.74 

0 

2.31 

4.94 

Total  C02d  (kgC02  kgH2-‘) 

8.99 

8.45 

7.71 

5.30 

4.99 

a  CH4  consumed  in  the  process  (Table  2)  or  as  hot  utility  (Table  1). 
b  CO2  emitted  from  the  process  or  as  hot  utility.  Dry  ice  not  included. 

c  CO2  emitted  (indirectly)  as  work  is  generated  in  a  power  plant  to  make  it  available  as  a  utility.  Estimated  using  a  conversion  factor  for  electricity  of 
0.43  kgC02  kW-1  h-1  [26]. 

d  Total  CO2  emission  =  plant  CO2  +  work  associated  CO2.  Dry  ice  not  included. 


deep  ocean  sediments  where  the  carbon  dioxide  is  chemically 
sequestered  via  the  formation  of  a  clathrate  (CO2 -hydrate) 
[25], 

Calculations  of  CO2  emissions  are  summarized  in  Table  4. 
Due  to  a  much  lower  consumption  of  CH4  (for  both  process 
and  hot  utility)  and  CO2  sequestration  in  the  form  of  dry 
ice,  the  heat  and  power  integrated  PE  based  process  offers 
a  great  opportunity  for  plant  CO2  emission  reduction  (up  to 
99%)  with  respect  to  the  PSA  based  process,  but  this  brings 
along  the  need  for  additional  work  that  has  the  potential  for 
an  associated  CO2  emission  if  fossil  fuels  are  used  in  the  gen¬ 
eration  of  such  work.  Assuming  a  conversion  factor  in  elec¬ 
tricity  generation  of  0.43  kg  CO2  kW-1  h-1  [26],  the  total 
CO2  reduction  is  35%  over  the  heat-power  integrated  PSA 
process. 


6.  Conclusions 

A  clean  hydrogen  (99.9999%)  and  dry  ice  production  pro¬ 
cess,  that  employs  methane  reforming  and  phase  equilibrium 
(PE)  based  separation,  has  been  proposed  and  simulated.  Heat 
and  power  integration  studies  have  been  carried  out  for  the 
proposed  process  in  order  to  determine  its  minimum  utility 
cost.  Using  this  optimum  value,  the  operating  cost  of  this 
alternative  process  is  then  calculated.  It  is  found  that  if  the 
produced  dry  ice  is  sold  for  as  low  as  2  cents  kg-dry-ice_  1  or 
if  an  equivalent  CO2  sequestration  credit  is  conceded,  the  PE 
based  process  results  in  lower  operating  cost  than  the  con¬ 
ventional  PSA  based  process. 

After  proper  disposal  of  the  dry  ice,  like  CO2  storage  in 
the  ocean,  the  PE  based  process  achieves  a  35-99%  reduc¬ 
tion  in  CO2  emissions  with  respect  to  the  PSA  based  process, 
depending  on  whether  or  not  fossil  fuels  are  used  for  elec¬ 
tricity  generation. 


Acknowledgement 

The  authors  gratefully  acknowledge  the  financial  sup¬ 
port  of  the  National  Science  Foundation  under  Grant  CTS 
0301931. 


Appendix  A.  Kinetic  models  for  reactions  ri,  r2  and 
r3  in  the  reformer 

Model  proposed  by  Xu  and  Froment  [15]  for  reactions  on 
a  Ni/MgAUCU  catalyst: 


k]  (  P('ii4 P\ho  P^Pco\ 

11  “  Dei?  ^  P2.5 

_  k.2  f  Pco  P[  [2p  _  I\:o2  \ 

'r2  “  Den2  V  Ph2  K2  J 

h  (  k\  \  i4  Pfho  P|°[25  Pco2  \ 

’I3  “  Den2  (  4.5  K3  J 


(A.l) 

(A. 2) 


(A. 3) 


where 


Den  =  1  +  kcoPco  +  &h2  Pr2  +  &ch4  Pch4 


■  Ph2o 


Ph2o\ 
Pu2  J 


(A.4) 


The  rate  coefficients  and  adsorption  constants  are  given 
in  Table  A.l  and  the  reaction  equilibrium  constants  in 
Table  A.2. 


Table  A.l 

Rate  coefficients  and  adsorption*  constants’  parameters  for  related  Arrhenius 


or  Van’t  Hoff* 

equations 

Rate 

Pre-exponential 

Units  of  pre-exponential 

Activation 

coefficient  or 

adsorption 

constant 

factor 

factor 

energy  or 
adsorption 
enthalpy 
(kJ  mol-1) 

h 

4.225  x  1015 

kmol  bar0'5  kgcat-1  h_1 

240.1 

1.955  x  106 

kmolkgCat_1  h_1  bar-1 

67.13 

h 

1.020  x  1015 

kmol  bar0'5  kgcat-1  h_1 

243.9 

0 

u 

8.23  x  10-5 

baq1 

-70.65* 

kcru* 

6.65  x  10~4 

bar-1 

-38.28* 

O 

X 

1.77  x  105 

dimensionless 

88.68* 

kn2* 

6.12  x  10-9 

bar-1 

-82.90* 
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Table  A.2 


Reaction  equilibrium  constants 


Equilibrium 

Function  of  T  (K) 

Units  of  equilibrium 

constant 

constant 

Ki 

exp(— 26,830/r  +  30. 114) 

bar2 

k2 

exp(4,400/r- 4.036) 

dimensionless 

Ki 

K\  x  Kj 

bar2 
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